A nonlinear symmetry breaking effect in shear cracks 
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Shear cracks propagation is a basic dynamical process that mediates interfacial failure. We develop 
a general weakly nonlinear elastic theory of shear cracks and show that these experience tensile-mode 
crack tip deformation, including possibly opening displacements, in agreement with Stephenson's 
prediction. We quantify this nonlinear symmetry breaking effect, under two-dimensional deforma- 
tion conditions, by an explicit inequality in terms of the first and second order elastic constants 
in the quasi-static regime and semi-analytic calculations in the fully dynamic regime. Our general 
results are applied to various materials. Finally, we discuss available works in the literature and 
note the potential relevance of elastic nonlinearities for frictional cracks. 



A fundamental process that governs the strength and 
failure of interfaces in a wide range of scientific and tech- 
nological problems is the propagation of shear cracks. For 
example, stress and energy release during earthquakes is 
mediated by the propagation of shear cracks along fault 
surfaces [l|. In generic situations, shear driven cracks 
(also known as mode II symmetry cracks [2.]) involve 
interfacial slip, but no crack tip opening displacements 
and/or tensile stresses, which are characteristic of mode 
I symmetry cracks [2]. However, the existence of such a 
symmetry breaking effect is of interest and potential im- 
portance. This becomes evident when the propagation of 
shear cracks is accompanied by frictional processes. Un- 
der such circumstances, crack opening, or even a reduc- 
tion in the confining normal stress at the interface, can 
significantly reduce the friction force and hence result in 
interface weakening. 

The most widely invoked mechanism for the coupling 
between interfacial slip and normal stress variations (and 
possibly crack opening) is material discontinuity, i.e. the 
shear crack is assumed to propagate along an interface 
separating two dissimilar materials [3]. An alternative 
physical mechanism, operative even in the absence of ma- 
terial discontinuity, is elastic nonlinearity 0, H| • 

In this Letter, we explore the latter mechanism by de- 
riving general weakly nonlinear asymptotic solutions for 
frictionless shear cracks in homogeneous isotropic mate- 
rials. The conditions for crack tip opening, as a func- 
tion of the first and second order elastic constants and 
crack propagation velocity, are quantified and the theory 
is applied to polymeric, glassy and metallic materials. Fi- 
nally, we discuss the possible relevance of our results to 
frictional shear cracks and to the super-shear transition 

We note that while the absence of pure shear crack so- 
lutions in nonlinear elastic fracture mechanics was recog- 
nized by Stephenson in 1982 [5], to the best of our knowl- 
edge no accurate and self-consistent benchmark weakly 
nonlinear solutions are presently available (see below). 

We start by considering the motion <p, which is as- 
sumed to be a continuous, diffcrentiable and invert- 
ible mapping between a reference configuration x and 



a deformed configuration x' such that x' = ip(x) = 
x + u(x), where u(x) is the displacement field. The 
Green-Lagrange (metric) strain tensor E is defined in 
terms of the deformation gradient tensor F = Vtp = I+H 
as E=^(F T F—I) '7J. Here I is the identity tensor and 
H = Vu is the displacement gradient tensor. 

The elastic energy per unit reference volume U(E) of 
isotropic and compressible materials can be expressed in 
terms of the scalar invariants of E, which are conve- 
niently written as trE, trE 2 and trE 3 . To third order 
in E, U(E) takes the form 



U(E) = - (trE) 2 + ^ trE 2 
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where {A, p} are the first order (Lame) constants and 
{«i}i=i-3 are second order elastic constants. The latter 
are important physical constants that are closely related 
to the leading anharmonic contributions to interatomic 
interaction potentials. These anharmonic contributions 
are known to be the origin of many basic physical phe- 
nomena and properties such as the Griineisen parame- 
ters, deviations from the Dulong-Petit law at high tem- 
peratures, the thermal expansion coefficient and the exis- 
tence of thermal resistance. Their implications for mate- 
rials failure remain largely unexplored (but see (HH)])) a 
situation that the present Letter aims to at least partially 
improve. 

For simplicity, in what follows we focus on plane strain 
conditions [3] for which Eq. ([1]) reduces to 
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where E is understood hereafter as a 2D tensor and 
I = ct\ + a 2 and m = a 3 + a 2 are two of the Murnaghan 
coefficients [7|. 

To develop a fracture theory based on this general ex- 
pansion we write down the momentum balance equation 
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where p is the reference mass density and s = OfU 20 
is the first Piola-Kirchhoff stress tensor, that is work- 
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conjugate to the deformation gradient F 0j. A friction- 
less crack is introduced by the usual traction-free bound- 
ary conditions 



±ir) = s yy (r, 9 = ±n) = , 
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where (r, 9) is a polar coordinate system that moves with 
the crack tip and is related to the reference frame by 
r = \/ (x — vt) 2 + y 2 and 9 — tan _1 [y/(x — vt)}. a; is the 
propagation direction, y is the perpendicular direction 
and v is the tip propagation velocity. In what follows we 
neglect crack acceleration effects, implying that all of the 
fields depend on x and t through the combination x~vt 
and therefore d t = —vd x . 

To proceed, we expand s order by order in powers of 
the displacement gradient H to obtain 



s(H) ~ A trel + 2pe + -A tr{H T H)I + fi H T H 
A treH + 2/i He + I {trefl + 2me 2 + 0(H 3 ) , 



(5) 



where e = \ {H + H T ) is the infinitesimal strain tensor. 
We then follow [1 - 15 , 17 1 and write the displacement field 
as 



u(r, 



(2), 
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where the superscripts denote different orders in H . Sub- 
stituting Eqs. ©-© into Eqs. ©-©, we obtain the 
equations and boundary conditions for and u^ 2 \ 

To linear order, we obtain the well-studied equations 
of isotropic Linear Elastic Fracture Mechanics (LEFM) 
0. To second order, we obtain 

M VV 2) + (A + M )V(V • m (2) ) + F[u {1) ] = pu (2) , (7) 

where the boundary conditions at 9 = ±ir read 



-(A + 2n) r- l d e u y V -\d r u^ - Sy[uW] = . (8) 



The functionals ^"[m*- 1 ^] and .Sl/u^] are easily obtained 
from Eq. ([5]). Asymptotic near-tip solutions of for 
mode I (tensile) symmetry were derived in [Tj 1 "i . with 
a focus on a particular constitutive law. 



Here we focus on global mode II (shear) symmetry 
loading. Under these conditions, the linear order solu- 
tion m' 1 ) possesses the following symmetry properties 

u£ ) (r,-6) = -u£'>(r,6) and ^'(r, -6) = u y 1) (r, 6). (9) 

The leading order near-tip asymptotic solution reads Q 
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pV2ttD(v) 

2a s a d r 1 d /2 cos (9 d /2) - (1 + a 2 )^ 1 / 2 cos (9 s /2) 
l - v2 / c d.s' tan #<ls = a dtS teav9, r dt 



Here a 2 s 

r^l-(vsin9/c d , s ) 2 , D{v) = 4a s a d - (1 + a 2 ) 2 ; K n 
is the mode II "stress intensity factor" and c d — 
y(A-j- 2/i)/p, c s = U \xj p are the dilatational and shear 
wave speeds, respectively |2|. 

Using Eqs. (JTOj) to calculate T[u^] and S[u (1) } in 
Eqs. CO)-©, we obtain 
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where g II {9\v) and k ti (v) are functions of stress di- 
mension that are too lengthy to report here (the super- 
script denotes that these quantities correspond to mode 



II loading conditions). J- (r,9) and «S y (r, ±7r), when 
substituted into Eqs. ©-©, can be respectively in- 
terpreted as an effective body force and a surface force 
in an effective linear elastic crack problem. Their sym- 
metry properties dictate the symmetry properties of the 
weakly nonlinear solution u^ 2 \ Equation (fTT|) indicates 
that «S 7/ (r, 7r) = S l:l (r,—ii) and an explicit calculation 
shows that 

g I x I (9;v)=gl I (-6;v) and gjfftv) = -<?^(-M • 

(12) 

These properties imply that the weakly nonlinear prob- 
lem is of pure mode I nature, even though the loading 
is of pure mode II nature |8l4lOj. Therefore, while it*- 1 ) 
possesses mode II symmetry, cf. Eq. ([9]), it*- 2 -* satisfies 

u£\r,-9) = uM(r,d) and uf\r, -9) = -u y 2 \r, 9), (13) 

which is the nonlinear symmetry breaking effect we aim 
at studying here. 

To quantify the effect we need to solve the problem 
posed in Eqs. ©, © and pTlK Since this is a pure 
mode I problem, we can follow [12h15| to obtain 
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T II (9;v) is an r-independent solution of Eq. ((TJ) that 
does not satisfy the boundary conditions in Eq. ((HJ. It 
can readily obtained semi-analytically using a Fourier se- 
ries method [12, [HI . The rest of the solution satisfies Eq. 
([7]) when the effective body force is omitted. It contains 
two parameters, A and B, that should be determined 
by the boundary conditions. However, the first equation 
in © is satisfied automatically and only the second one 
imposes a constraint of the form 



A = 



2fiBa s — (A + 2ui)d e T I y I (n- v) ~ n H (v) 
A- (A + 2/z)a 2 , 



(15) 



To complete the solution, we need an additional con- 
dition that allows the calculation of both A and B. To 
see where this additional condition emerges from, substi- 
tute the solution in Eq. (TUT) in s of Eq. (|5|) and write 
s = s^+s^ 2 \ where the superscripts denote different 
orders of Ku. is characterized by a 1/r spatial sin- 
gularity that was shown to give rise to a (spurious) force 

(2) 

in the crack parallel direction, f x , a force that cannot be 
balanced by material inertia which vanishes to this order 
HI HIGH- Therefore, we need to impose the condition 
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ensuring that our solution is consistent with the conser- 
vation of momentum, which is otherwise violated. Equa- 
tion (|16[) . which is a linear relation between A and B, 
together with Eq. (fTS"j) . allows us to fully determine the 
solution in Eq. (fM)) . 

A natural and physically transparent way to quantify 
the symmetry breaking effect is through the dimension- 
less crack tip displacement defined as 
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Crack tip opening emerges when A > 0. The com- 
plementary situation, A < 0, obviously does not imply 
that the crack's faces penetrate each other, but rather 
that the traction-free boundary conditions are no longer 
valid and other physical effects should be taken into ac- 
count. Equation (|14p immediately tells us that A is 
r-independent. Hence, we are mainly interested in the 
dependence of A on the first and second order elastic 
constants and on the crack velocity, i.e. A(A, p,,l,m,v). 
To study this function we first focus on the quasi-static 
limit v — > 0, and later examine the behavior in the fully 
dynamic regime. 



The problem for v = can be solved analytically and 
leads to 



A(t> = 0) = 

(l-2i/) 2 {l-2vfl i/(l - v){l - 2v) m 



(18) 



where v = A/[2(A + /i)] is Poisson's ratio. We first note 
that A depends only on three dimensionless material pa- 
rameters, v, I /n and m/fi. In addition, A contains a 
common factor inversely proportional to the bulk modu- 
lus (1 — 2v)~ 1 and may possibly feature an interest- 
ing behavior as the incompressibility limit is approached, 
v—> 1/2. We are mainly interested in the sign of A, and 
if positive, in its magnitude. Therefore, Eq. (|18j) can 
in fact serve as an explicit inequality that v, l//i and 
ml \x should satisfy in order to observe crack tip opening, 
A>0. 

In order to better understand this inequality, let us 
consider the different terms in Eq. (fT8")l . The first term 
is independent of the second order elastic constants I and 
m, and hence is associated with a geometric nonlinear- 
ity. This term exists even in the absence of constitutive 
nonlinearities, / = m = 0, and emerges because E is non- 
linear in H , see above. Most importantly, it is positive 
independently of v. Therefore, we conclude that geo- 
metric nonlinearities universally tend to induce opening 
of shear cracks in the quasi-static regime, v = 0. The 
last two terms on the right-hand-side of Eq. (fT8"|) de- 
pend on constitutive nonlinearities quantified by / and 
to. Since I and to can be both positive and negative 
(typically Z,m 3> /i), we cannot say something definite 
about these contributions and whether the inequality is 
satisfied or not depends on the specific values for each 
material. 

It would be interesting to analyze the dependence of 
A on the crack propagation velocity v. To extend Eq. 
(fT8| to the dynamic regime we should introduce another 
dimensionless parameter, v/c s , and write 

I TYl 

k{v) = C {v,v/c s )+C l {v,v/c s )~+C rn {v,v/c s )— , (19) 

where Co, C; and C m are dimensionless functions that 
reduce to their counterparts in Eq. ([TBI when v — > 0. 
A(v) was calculated semi-analytically following the pro- 



cedure described in |12Ml5| . In Fig. [T]we show A(i>) for 
a polymer (Polystyrene with v — 0.34, 1//J, = —13.7 and 
m/fi = — 9.6 [19jn. a glass (Pyrex with i/ = 0.17, Z//i = 0.5 
and to//i = 3.4 |l9|) and a metal (Copper with j/ = 0.35, 
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FIG. 1. A as a function of v/c s for three materials. 

l/li = 11.8 and m//i = -8.1 [20]). For Polystyrene, A(u) 
remains nearly v-independent until it changes sign at 
high v. For Pyrex, A(u) increases monotonically and 
quite significantly at high v. For Copper, A(v) increases 
up to a high v where it abruptly decreases and eventually 
becomes negative. The interplay between materials pa- 
rameters and elastodynamic effects gives rise to this rich 
behavior. 
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FIG. 2. A linear-log plot of the coefficient functions 
Co,Ci, Cm, defined in Eq. (fT9|). vs. v for n = 0.8c s . 

The extreme elastodynamic limit, v^c s , is of partic- 
ular interest. This regime might be relevant for phe- 
nomena such as ordinary earthquakes and the super- 
shear transition [6]. To better understand what deter- 
mines the sign of A at high v, we present in Fig. [2] 
Co, C, C m as a function of v for v — 0.8c s . Several 
qualitative conclusions can be drawn. All C"s are pos- 
itive and satisfy Co ^> C m > Cj, where C m ^> Ci for 
large values of v (y > 0.3). Therefore, I typically plays 
a minor role in determining the sign of A(w = 0.8c s ). In 
addition, positive values of m and I immediately ensure 
A(v = 0.8c s ) > 0, which is the case of Pyrex, cf. Fig. 
[TJ Finally, since Co is an order of magnitude (or more) 



larger than C m , a large and negative m is needed in order 
to have A(v = 0.8c s ) <0. This is the case of Polystyrene, 
cf. Fig. [Q 

To put our results in the context of existing litera- 
ture, we note that there were very few attempts to derive 
specific weakly nonlinear crack tip solutions |8h10|. In- 
variably, all of these attempts assumed that u^(r, 9) is 
r- independent. Equations (j!4[) clearly demonstrate that 

( 2 ) 

this is not the case as Mi (r, 9) contains a contribution 
proportional to logr, which was experimentally verified 
for mode I fracture [12l4l4| - More generally, Eqs. (fT4|) 
show that u^ 2 \r,9) cannot be separated into a prod- 
uct of an r-dependent function and a ^-dependent func- 
tion. In essence, the assumption of an r-independent 
u^ 2 \r,9) amounts to implicitly adopting the condition 
A + Ba s = 0, which annihilates the logr term in Eqs. 
(JT4J) . With this conditions at hand, indeed a solution 
of Eqs. (0, ([8]) and (fTTj) can be found. Unfortunately, 
such a solution does not satisfy Eq. (flo) . Therefore, the 
solutions discussed in produce an unphysical force 

in the crack propagation direction and hence violate the 
conservation of momentum. 

We have verified the latter statement in relation to 
the solution presented in Q, which was derived for a 
nonlinear Hookean material (I — m = in Eq. ([2])) in 
the quasi-static limit (v— >0). We recovered this solution 
using A+B = (recall that a s = 1 for v = 0) instead of Eq. 
(IT51) . However, as explained above, this solution violates 
Eq. (1T61) . i.e. it produces a finite fx ■ In addition, 
the solution in Q predicts crack closure, A < 0, which 
is inconsistent with Eqs. (|T8"1) that predicts A > for 
I = m = 0. 

The scale of the opening displacement A and the dis- 
tance from the crack tip in which it takes place are in- 
herited from the lengthscale Kjjf/j 2 . For example, for 
a material with /i = 1 GPa and Kjj = 5 MPa-^/m (which 
may be relevant for polymers) we obtain a lengthscale of 
~ 25/im. The value of Kjj (the dynamic mode II fracture 
toughness) is a material parameter that is, unfortunately, 
not well documented for many materials under dynamic 
conditions. All of our results can be recast in a dimen- 
sional form once Kjj is known. 

The symmetry breaking effect we have quantitatively 
addressed in this work may have physical consequences 
for several important problems. As was already men- 
tioned in the introduction, this effect may influence the 
strength of frictional interfaces. Since the local frictional 
resistance along interfacial contact depends linearly on 
the compressive normal stress, crack tip opening or even 
a reduction in this stress, may locally weaken the inter- 
face. Elastic nonlinearities, which are rather generally 
overlooked in frictional analysis, may provide a mecha- 
nism for such a local weakening near the crack tip. A 
quantitative analysis of this possible effect goes beyond 
the present work, entailing the introduction of a friction 
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law acting on the crack faces and the inclusion of a com- 
pressive normal stress. 



Another relevant problem is the super-shear transi- 
tion beyond which a shear crack propagates at a velocity 
larger than the shear wave speed [2l|. The transition 
from a sub-shear state to a super-shear state is believed 
to be mediated by the nucleation of a "daughter" crack 
in front of the "mother" crack [(|. Above, we quantified 
the symmetry breaking effect as manifested at the crack 
faces, i.e. 8 = ±7T. However, tensile components of the 
displacement and stress fields emerge also at a continu- 
ous range of other angles, including ahead of the crack 
tip, potentially facilitating the super-shear transition. In- 
deed, the numerical simulations of [gj], which demon- 
strated the super-shear transition, also feature opening 
displacement that might be interpreted as the symmetry 
breaking effect we addressed [loj . 



Wc hope that the present work, beyond providing gen- 
eral and self-consistent weakly nonlinear shear crack solu- 
tions, will serve as an impetus for more systematic studies 
of the roles played by elastic nonlinearities in the dynam- 
ics of shear cracks in various of problems. 
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